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Abstract. - We investigate the effects of magnetic impurities in a superconducting state with s± 
pairing symmetry. Within a two-band model, we find that the intra-band magnetic scattering 
serves as a pair breaker while the inter-band magnetic scattering preserves pairing and hardly 
affects transition temperature in the Born limit. We also show that the same physics can persist 
beyond the weak scattering region. Our results coincide with recent experimental measurements 
in iron-based superconductors and thus provides an indirect evidence of the possible s± pairing 
symmetry in these materials. 
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O . The iron-based superconductors have attracted much 
i ^ a ttention since the compound LaFeAsO\- x F x was found 
superconducting with T c = 26K [1]. Accompanied with 
C^J the increasing transition temperature over 40K [2] and fi- 
^nally up to 55K [3], these new superconductors are con- 



oo 



sidercd as the second family of high-tempcraturc super- 
QQqonductors after the cuprates. Though much efforts have 
(f~)been made and great progress has been achieved experi- 
ly-^ ipentally, the superconducting mechanism and the pairing 
(**••*) symmetry are still unclear because of the complexity of 
multi-orbital nature and possible strong electron-electron 
^^correlation in these materials. The angle-resolved pho- 
£>toemission spectroscopy (ARPES) [4], the Andreev reflec- 
KJi ^ on [5] an d the penetration depth [6] experiments directly 
rS sjhowed fully gaped superconductivity while the unclear- 
ed magnetic-resonance (NMR) indicates a strong deviation 
" from that of single band s-wave superconductors [7] . Si- 
multaneously, the s± pairing symmetry was proposed and 
widely discussed by many authors theoretically [8]. With 
this kind of pairing symmetry, fully opened gap located in 
the hole pocket around the T point and electron pocket 
around the M point have opposite signs. The s± pairing 
symmetry is consistent with the ARPES results and also 
coincides with the NMR results providing the presence of 
strong impurities [9] and thus becomes a very promising 
candidate to account for the main physics of the iron-based 
superconductors . 

To distinguish the pairing symmetry experimentally, 
we note that the effects of impurities in superconducting 
states can be very different for different pairing symme- 



tries. Very recently, superconductivity in both 1111 and 
122 systems has been induced by doping magnetic ele- 
ments Co [10,11]. These experiments share some common 
features. First, the superconductivity in these materials 
shows high tolerance with the disorder induced by the Co- 
doping. Second, the suppression of the transition temper- 
ature is not so significant but still much stronger than the 
F-doped materials. The second point can be easily under- 
stood because the F-doping happens in the LaO layers 
while the Co-doping happens in the FeAs layers which 
are crucial to the superconductivity. However, a later ex- 
periment [12] on Zn-doped LaFeAsO shows that the su- 
perconductivity is almost unperturbed by the Zn-doping, 
though it happened in the FeAs layers. This strongly indi- 
cates that the non-magnetic impurities are unlikely to af- 
fect the superconductivity, in accordance with s-wave pair- 
ing symmetry. Meanwhile, magnetic impurities can de- 
press the transition temperature but not so significant as 
in the conventional s-wave superconductors. A. E. Karkin 
et al used fast neutron irradiation to induce disorder in 
FaFe/lsOo.9-F).i and found the depression of transition 
temperature can be qualitatively described by Abrikosov- 
Gorkov (AG) theory with magnetic impurities [13]. How- 
ever, such depression is also much slower than that pre- 
dicted by the AG theory (see Fig. 10 in [13]). 

To solve this puzzle, we propose a model in this let- 
ter to describe the behavior of the magnetic impurities in 
an s±-wave superconductor and show how the inter-band 
magnetic impurity scattering processes may preserve the 
pairs and effectively weaken the reduction of transition 
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temperature. Our theoretical results give a reasonable ex- 
planation of the experimental data. 

We start with a model consisting of two perfect nested 
bands, i.e., an electron Fermi pocket and a hole Fermi 
pocket. The details of the band structure is neglected 
and the superconducting order parameters located in each 
band have same magnitude but reversed signs as in Refs. 
[14, 15]. For convenience, we introduce the Nambu vector 
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here c k(T and d k(T are the annihilation operators in the 
electron band and the hole band and the band energy 
e Cjk = — £d, k = e k , respectively. Here we treat the mag- 
netic impurities as localized spins in the classical limit 
and the quantum (Kondo) effect of impurities is not un- 
der our consideration. In this limit, the magnetic impu- 
rity is equivalent to the local magnetic field. Then the 
interaction matrix due to magnetic impurity scattering is 
assumed to be 
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Here a and S denotes the spin operator of the electrons 
and magnetic impurities, respectively. Both the intra- 
band (Ji) and the inter-band {J 2) exchange coupling con- 
stants are positive and isotropic in our consideration. 

Following the AG theory [16], the renormalizcd two- 
band BCS Green's function with randomly distributed im- 
purities of finite concentration rii mp reads: 
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where Go(k, lu) is the Green's function without the im- 
purity; XX W ) i s the self energy; uj and A is the renor- 
malizcd frequency and superconducting order parameter, 
respectively. The parameter x is the inter-band scattering 
induced contribution to the self-energy which can be de- 
termined self-consistently. In the Born approximation the 
self-energy can be written as [17] 
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where (...) j means averaging the impurity position. Now 
the renormalizcd Green's function can be given self- 



consistently with 
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and the spin-flip scattering times: 
1 (ff + Jl), 
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with Np the density of states at the Fermi surface. 
With the above equations we obtain 

— = u[l - 

T S A Jl + U 2 - (1)3 



with u = ~ and I = ix. The effective pair-breaking pa- 
rameter is a = = (i o^- + tt- + tM4-. 

Now it is clear that the Green's function contains two 
terms: one is the conventional term as in the single- 
band s-wave superconductors with the renormalizcd fre- 
quency and superconducting order parameter and the 
other one is the contribution induced by the inter-band 
magnetic impurity scattering processes. If the inter-band 
scattering term x in the renormalizcd Green's function 
is neglected, the effective pair-breaking parameter reads 
a = cto = (jij- + 2^)s ^ "^l' quite similar to that in 
the conventional AG theory. However, a non-zero J2 will 
change the situation significantly as we shall show below. 

Making iuj — > ui, the density of states (DOS) is given 
by: 

N(uj) = --Im [ -^LG*(k,u) 

IT J (27TP 



Nplm( 



V A 2 - uj 2 + x 2 



)■ 



Numerical results of the DOS for given ao and A = J 2 /Ji 
are shown in Fig.l. The densities of states for both the 
conventional two-band s-wave case and the s±-wave case 
are calculated. The numerical results clearly show that 
in the conventional s-wave case (left part of Fig.l), the 
inter-band scattering also depresses the superconducting 
gap. However, accompanied with the increasing of J2, the 
superconducting gap is growing larger in s±-wave case. 
This strongly indicates that J\ and J2 have opposite effects 
on the gap of s±-wave superconductors. It seems that the 
inter-band magnetic impurity scattering played as a pair 
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Fig. 1: (Color online) The density of states for ao = 0.2 and 
A = J'i/Ji- The left part is for the s-wave pairing symmetry 
and the right one is for s±-wave pairing symmetry. 



repairer in our system. To make that clearer we calculated 
the superconducting gap in finite temperature. Here the 
superconducting gap A in equation (3) should be replaced 
by A(T), and A(T) is determined by 



A(T) = V sc N F nT^2 



A 



(5) 



where V sc is the coupling constant and oj m = (2to+1)7tT. 

Solution of (2)-(5) gives the finite temperature gap 
A(T). The intra- and inter-band magnetic impurity ef- 
fect to A(T) are shown in Fig. 2. In Fig. 2(a) only intra- 
band impurity scattering exists and — = ~ oc J\ . When 
we increase the intra-band scattering, the finite tempera- 
ture gap and the transition temperature become smaller. 
This case is in accordance with conventional s-wave su- 
perconductors [18]. In Fig. 2(b) things will be reversed if 
we settle down the intra-band scattering and increase the 
inter-band one. A(T) and Tc becomes larger with J2 in- 
creasing. 

Such an effect also reflects in the transition temperature 



7 Tc 
In 

T 

1 r.r> 



m=0 



1 



1 



(™ + 5K/i-(£) 2 + 



rs 2ttT c 



where T c and T cp are the transition temperatures with and 
without magnetic impurities, respectively. The numerical 
results of the transition temperature vs. ao is depicted 
Fig. 3. Once the A = J2/J1 increasing, the depression of 
the transition temperature is effectively speeded up in the 
s-wave superconductor but slowed in the s±-wave super- 
conductor. In experiments, Co-doping introduces extra 
carriers and modifies the crystal structure [10]. As those 
changes may be crucial to the superconductivity, it is quite 
difficult to examine whether other issues may play a role 
in the robustness of superconductivity. However, a very 
recent neutron irradiation experiment on LaO0.9F0.1Fe As 
provides a chance to check our theory since in this experi- 
ment both the crystal structure and carrier density of the 
sample are almost unchanged. The comparison is shown 
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Fig. 2: (Color online) The superconducting gap in finite tem- 
perature with different intra- or inter-band impurity scatter- 
ing amplitudes. The solid black line in both (a) and (b) is 
the two bands BCS case without impurities. A p (0) is the 
superconducting gap without impurities in zero temperature. 
(a)J2 = 0. (b) The solid black line: rii mv = J\ = 0; others: 



,N F S(S + 1) =0.7, Ji = 1. 



in Fig. 3. The experimental data is fitted quite well with 
our theory quantitatively. In the real iron-based materi- 
als, with the increasing of the Co-doping, the hole pocket 
becomes smaller and the electron pocket becomes larger 
because of the shift of the chemical potential. This will 
weaken the inter-band magnetic scattering and preserving 
of T c . 

If A > 1, T c goes down to a finite value rather than 
zero when ao is very large. This indicates that in the 
Born approximation the superconductivity can never be de- 
stroyed by magnetic impurities if the inter-band scattering 
is stronger than the intra-band one in the s± -wave su- 
perconductors. To clarify this issue further, we consider a 
single classical spin in an s±-wave superconductor, which 
can be treated with more accuracy. The T-matrix in this 
case is given by 

T = V[l-G°(0,u)V}- 1 . 

The energy of the bound state induced by the single mag- 
netic impurity in an s± superconductor is determined by 
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Fig. 3: (Color online) The reduced transition temperature for 
the conventional s-wave (above) and s±-wave (below). When 
A = the conventional AG behavior is regained. Inset: fitting 
of the experimental data (filled squares) from Ref. [13] with 
A = 0.93 



the pole of the T-matrix 
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while the position of the bound state for the conventional 



s-wave superconductor is ^ = ±{ 
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If J 2 = 0, i.e., only the intra-band impurity scattering ex- 
ists [19], the bound state energy (BSE) is 



"0 

A 
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1+ai 



which falls into the gap [20] and we recover the Yu-Shiba- 
Rusinov solution [21-23]. If only the inter-band impurity 
scattering exists, The BSE is ^ — > ±1, i.e., locates at the 
gap edge. Generally, the BSE falls into the gap. Note for 
the conventional two-band s-wave superconductors, the 
bound state splits into two branches in each band due 
to the inter-band magnetic scattering. However, for the 
s± case, there is only one bound state and increasing J\ 
pushes the BSE to the Fermi energy side while increasing 
J 2 pushes the BSE to the gap edge side. With a finite 



impurity concentration, an impurity band [22] will be ex- 
panded around the position of the BSE. The supercon- 
ducting gap can be suppressed (enlarged) by increasing 
J\ (^2)- This analysis supports our Born approximation 
result. 
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Fig. 4: (Color online) Left: The position of the bound state 
energy vs. the inter-band coupling constant J2 in conventional 
s and s± superconductors. The two branches of the bound 
states in s case are labeled as BSE1 and BSE2 and Ji = 0.75. 
Right: Change of the position of the bound state energy with 
Qi and A = J2/J1 

When A > 1, from Eq.(3) we can see that 

1 — OLl + a 2 
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which means in this case the BSE has always a positive 
value and can never reach the Fermi level. Such fact pre- 
serves the finite superconducting gap even with a finite 
impurity concentration of impurities and explains why T c 
goes to a finite value rather than zero with the increasing 
of ao in Fig. 3. The position of the BSE with increasing a\ 
and different A is shown in the right part of Fig. 4. Con- 
sidering the impurity band expanded around the position 
of BSE [22], the region where always exists a finite gap 
should be in A > 1. The discussion on the single impurity 
problem makes our conclusion to a broader region beyond 
the weak scattering limit. 

It is well known that the magnetic impurity is a pair- 
breaker in the conventional spin singlet superconductors 
because it breaks the time-reversal symmetry and non- 
magnetic impurity can break the pairs in e?-wavc super- 
conductors because it smears the anisotropy of the order 
parameter. If we neglect J\ and look back to the effect of 
the inter-band impurity scattering, we find that both of 
the time-reversal breaking and sign reversal of the order 
parameters exist but canceled each other in the s±-wave 
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superconductors. In fact, due to the sign reversal of the 
order parameters, the wave functions in different pockets 
have a spin up-down exchange. Therefore, the inter-band 
spin-flip scattering process caused by magnetic impurity is 
the same as the intra-band nonmagnetic scattering process 
which preserves the spins and does not break the pairs. 
This is in accordance with the "reversed AG theory" in 
the inter-band channel which is proposed and noticed be- 
fore [24] . Though the inter-band magnetic scattering has 
no effect on pair breaking but J\ can still suppress the su- 
perconducting gap and T c . In our system which is based 
on one electron band and one hole band with the s±-wave 
pairing symmetry, the emerge of the inter-band scattering 
via impurities is inevitable and makes the inter-band mag- 
netic scattering behaves like a "pair repairer" and weakens 
the depression of the superconductivity. 

Based on our theoretical prediction, one may distin- 
guish the pairing symmetry in experiments by substitut- 
ing magnetic ions out of but coupled to the FeAs lay- 
ers and the T c depression must behaves quite differently 
in s-wave and s±-wave superconductors. A very recent 
experiment [25] showed that BaFe\- x Co x As2 is always 
more robust and have larger superconducting region than 
BaFe\- x Ni x As2, though Ni doping can result in smaller 
c-axis and bring more extra electrons into the supercon- 
ducting layer. In another recent experiment [26], the com- 
pound Ba\- x K x FeiAs-i shows a very large superconduct- 
ing area with doping (from x=0.f to x=1.0). These experi- 
ments also indicate the importance of the inter-band mag- 
netic scattering. However, the in-plane doped magnetic 
impurities may affect other factors relevant to the super- 
conductivity which make the situation unclear. One can 
also distinguish the s± and conventional s pairing symme- 
try in two band superconductors by detecting the bound 
state energy. There is one bound state in s± case but 
two in s case. Besides, with the increasing of J2, the BSE 
will move to the gap edge side in s± case. In conven- 
tional s case, BSE1 will move to the the Fermi edge side 
and quickly a quantum phase transition happens [27], but 
BSE2 first moves to the gap edge when J2 < J\ and then 
goes to the Fermi energy side when J2 > J\. Such be- 
havior is shown in the left part of Fig. 4. The position of 
BSE near the magnetic impurity can be detected from the 
tunneling spectra with a low-temperature scanning tun- 
neling microscope (STM) which has been used to detect 
the Yu-Shiba-Rusinov bound state successfully [28]. 

In summary, the magnetic impurity effect in the two- 
band superconductors with a perfect nesting effect is stud- 
ied for both the conventional s-wave pairing and the s±- 
wave pairing. Under the condition of the same magnitude 
of order parameters in each band, it is found that the de- 
pression behaviors to the superconductivity with these two 
different pairing symmetries are quite different due to the 
existence of the inter-band impurity scattering, which is 
almost inevitable in the multi-band systems. Our theory 
surprisingly coincides with the neutron irradiation experi- 
ment quantitatively and we believe it can also explain why 



the iron-based superconductors show high tolerance with 
magnetic impurities in many recent experiments. This 
provides an indirect method to detect the pairing symme- 
try of the FeAs superconductors in experiments. 
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